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spirit using the scattering wave functions. Based on this formulation, we derive some rigorous results 
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sandwiched symmetrically between two leads. These symmetry property can provide accuracy tests 
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I. INTRODUCTION 

During the past two decades, there is a fundamen- 
tal progress in the understanding of transport proper- 
ties in mesoscopic systems Quantum coherence in 
these systems plays a much more important role than 
in macroscopic dissipative systems. Accordingly, tradi- 
tional approaches based on quasi-classical picture (e.g., 
Boltzmann approach) of carriers ceased to work in these 
systems. Fortunately, a well known new approach was 
built up to take into account quantum coherence ade- 
quately, i.e., Landuer-Buttiker theory(LBT)0. In this 
theory, a mesoscopic system under consideration is con- 
nected to electric contacts through some ideal leads and 
the transport properties are determined by carrier scat- 
tering probability between these leads. In practice, the 
central part of LBT is the quantum mechanical scattering 
problem for a particle incident from each lead. However, 
under most circumstances, only the terminal quantities 
are actually concerned (most real experiments probe only 
terminal currents instead of current density, or any other 
local quantities in the system). Then, all one need is 
the scattering matrix between the leads and the terminal 
conductances can be simply cast into a beautiful formula 
expressed by Green's functions [lj. This formula is fre- 
quently used in literature. Up to our knowledge, no one 
has bothered himself previously to solve the whole scat- 
tering wave function explicitly. However, as is easy to 
see, the scattering wave functions are needed, if we want 
to calculate local properties in the Landauer-Buttiker 
Scenario. 

Recently, there is a hot topic in the field of spintron- 
ics, called spin hall effect(SHE)@, H ||. SHE refers to 
the phenomena that when a longitudinal electric field 
(or charge current) is applied in a semiconductor strip, 
transverse spin current and/or ensuing spin polarization 
near transverse boundaries can be induced. Such an ef- 



fect opens a possible new way to manipulate spin de- 
gree of freedom by all-electrical means, which is a big 
goal to the frontier research community in semiconduc- 
tor industry0,Q. Some years ago, two theoretical works 
predicted such an effect [j, H as an intrinsic property of 
spin-orbit coupled semiconductor band structure. Since 
then, a lot of theoretical efforts are made to clarify some 
fundamental controversies in this field0JHJSG3[llll3 
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are now also two central experimental results reporting 
the local measurement of opposite spin polarization near 
two transverse boundaries |27l l28j . which have attracted 
many theoretical interpretations. Despite many impor- 
tant progresses and consensus made, there sill remain 
some basic difficulties pif. 

Different theoretical approaches for spin dependent 
transport arc employed in this field. Among them, the 
Landauer-Buttiker formalism is suitable to address the 
ballistic transport regime and has the merit of fully tak- 
ing into account the phase information. This approach 
is in some sense mutually complementary to Boltzmann 
approaches based on the quasi-classical wave packet pic- 
ture. Based on this approach, some important numer- 
ical results has been achieved in the topic of intrinsic 
SHE[3 13 13 13 13 ES El El However, it is 
worthy to point out, that most previous works based 
on the Landauer-Buttiker's picture does not solve the 
scattering wave functions explicitly 0, 0, 0, 0, E3- 
The normal green's function formula[lj compute terminal 
conductances but not local quantities such as spin den- 
sity. On the other hand, the so called Landauer-Keldysh 
formalism [T2, fl3t |20| used by some authors which can 
compute local quantities is not manifestly single particle 
description and lacks intuitional simplicity. So, in or- 
der to calculate the non-equilibrium local quantities in 
a transparent manner, it is appealing to solve the sin- 
gle particle scattering wave function in the Landauer- 
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Buttiker geometry explicitly. 

In this paper, we will set up the linear equations to 
solve the scattering wave function in the multi-terminal 
Landauer geometry explicitly. With single particle wave 
function at hand, we can calculate any non-equilibrium 
local quantities(e.g.,spin density), by taking proper en- 
semble average in the Landuer-Buttiker spirit. Further- 
more, we will derive some rigorous results in the scat- 
tering wave function description. These include: (^de- 
duction of the green's function formula for terminal spin 
current widely used in literature, (2) rigorous proof of 
some important properties for equilibrium state in the 
Landauer's structure, (3) revelation of some new symme- 
try relations in a typical two terminal structure. These 
symmetry property can provide accuracy tests for exper- 
imental measurements and numerical calculations. 



II. DESCRIPTION OF THE METHOD TO 
SOLVE SCATTERING WAVE FUNCTION 
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nearest pair across the lead/system interface. H sys — 
Hq + H so is the Hamiltonian for the SO coupling system 
in which H = -tJ2 <Rz ,R J>a ( C R i a C R^+ h - c -) is the free 
band Hamiltonian, where Cr^o is annihilation operator in 
the SO system, and H so denotes SO coupling term. For 
a typical Rashba model, we write: 



ttR _ 
He„ — 



o-y* Ri+x ) + h.c] 



(2) 

Where t R is the Rashba coupling strength and 1 ^R i = 
(Cfl;,T' C-Ri.i) i s the vector representation of spinor anni- 
hilation operators. 

In the following, we will adopt local coordinate frame. 
Generally we use double coordinate index (x p , y p ), x p — 
1, ...co and y p = 1, N p for sites in lead p, and (x s , y s ) 
for sites in 2DEG system.We will also use single coordi- 
nate y or y s to represent one site on the leads or in the 
system when there is no risk of confusion. Furthermore, 
we denote n py as the boundary site in the system near 
to first row site (x p = 1, y p ) in lead p. 

Now, suppose there is an incident wave 

e mXp Xma(yp) m l ca< i P ( we assume the m'th 
transverse mode in lead p is conducting for incident 
energy E), where fc^ is the longitudinal wave vector, 
Xm is the wave function of m'th transverse mode and 
a is the spin index of the incident carrier. Hereafter 
we choose z axis (normal to the 2DEG plane) as the 
quantization axis for spin. Generally, the wave function 
in lead p' can be written as: 



(3) 



FIG. 1: schematic geometry of multi-terminal structure 



where 4^,'m'a-' i s the scattering amplitude from in-going 
mode (ma) in lead p to out-going mode (mV) in lead 
p' . The longitudinal wave vectors is determined by: 



We will consider a general mesoscopic structure in 
which a spin-orbit coupling system is attached to several 
ideal leads, as shown in Fig.l. The p'th lead is extend- 
ing to reservoir with chemical potential \i p at infinity. In 
a discrete representation, both the SO system and leads 
are modelled by a nearest neighbor tight binding(TB) 
Hamiltonian. We assume a tunable coupling between 
p'th lead and SO system by hopping interaction t ps . The 
Hamiltonian for the total system reads: 



H = H, 



leads 



H. 



sys 



Hi. 



(1) 



Here, Hi 



leads 



<Pi,Pj> 



I'itjCpjO 



+ h.c.) 



is 



Pj cr 



the TB Hamiltonian for the half infinite leads and C 
denotes annihilation operator for j 'th site with spin index 

a on p'th lead. H int = - J2 P tps S„(Cp„<jCfl„o- + h - c -) 
describes coupling between SO system and leads, where 
tp' s is the tunable nearest-neighbor hopping parameter 
between the lead p 1 and the SO system. (p n ,R n ) is the 



k P ' 



E - F P 



E-el 



< 2U 



(4) 



In the Eq.Q, the first line describes real fc^, and 
the corresponding mode is conducting. The last two 
lines have imaginary k v m , value and the corresponding 
mode is non-conducting, or, evanescent. The evanes- 
cent mode describes local components of scattering wave 
function in leads which decay exponentially away from 
the lead/system interface. It contributes zero to termi- 
nal current. The incident energy E lies out of the band 
continuum of such modes. 
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Now let's use a column vector ip s of length 27V (N is 
the number of lattice sites in the SO system) to repre- 
sent the wave function in the SO system ^T" 7 (x s , y s ) 
as a whole. Furthermore we arrange the scattering am- 
plitudes ^p'm'a' m t° another column vector <j> of length 
2M ( M = J2 P N p is the total number of transverse sites 
on all leads). Of course, these two vectors are the cen- 
tral quantities in our scattering problem and we want 
to solve them from the Schrodinger equations. As dif- 
ferent from the continuum model, the Schrodinger equa- 
tions equations has a lattice form( in real space, there 
is one separate equation centered on each site and spin 
index ) in our case and the normal boundary conditions 



has different appearance. Form Eq.© is the general lin- 
ear combination of scattering eigen-states in the leads, 
thus for all sites on the leads except the first row sites, 
the Schrodinger equation is satisfied automatically. The 
Schrodinger equation for first row sites, however, has a 
different form due to its coupling with the system. On 
the other hand, all sites in the system has normal form 
of Schrodinger equation for a closed system except the 
boundary sites which has coupling to the first row sites 
in the leads. These two connecting conditions constitute 
the boundary condition in our problem. To be explicit, we 
write down Schrodinger equations centered on each site 
in the system and on first row of all leads separately as: 



y'Ja" p',y 

y" Vs 



(5) 



Note on the above equation we have used y' s , y" to denote 
a site in the system and y' p , , y'L to denote a site in lead 
p' . However, at the same time we also used double index 
0->y' p ') etc. for first row sites in lead p' for clarity. n p / y i 
denotes the boundary site in the system near to first row 
site of lead p'. Substituting Eq.@ into Eq.JSJ) 

and after some routine algebra, wc arrive to the following 
matrix equations for the unknown vectors ip s and </>: 

AV> S = b + B(j) (6a) 



Ccj) = d + DV> S (6b) 

Here A, B,C,D are matrices with dimensions 2N x 
2JV, 2JV x 2M, 2M x 2M, 2M x 2N respectively. They 
are determined by the Hamiltonian and the geometric 
structure of the entire system, b, d are constant vectors 
determined by incident wave. These equations reflect 
the mutual influence of the leads and system through 
interface hopping. With some algebra, we can write down 
the matrix elements explicitly as following: 

A = E — H sys 

B(n p » y a" ,p'm'a') = - V'.p'^-'V VsXm'(%>") e m ' 

T>(p'm'a',np" v a") = -S p » yd<j» , CT 'VsXm'(%>") 

C(p'mV \p"m"a") = -<V>"V', CT '£ m "m' V 

b(n p/y a') = -5 pp >5<j<j>tpsX p m (yp') e ~ tkVm 

d(p'mV) = 5pp'5 mm '5aa't p (7) 



The indices for lead, transverse mode, site and spin 
can take all possible values. All other matrix elements 
not included in Eq.JTJ) are zero. Note for some corner 
points in the SO system, there may be nearest neighbor 
points in different leads. These corner sites have multiple 
notation in the above scheme. However, it causes no 
problem because in real programming, of course, we will 
assign certain index to one point. 

Note matrix C is diagonal and can be simply inverted. 
Defining matrix E fl = BC _1 D, it is easy to see that its 
elements are: 

^K lp ,a' )%2p ,a') = -5:f^X^(2/i)x6(2/ 2 )e<'' 

m' ^ 

(8) 

and all other elements are zero. From Eq.© we have: 
^ = (A-SV(b + BC- 1 d) = G R g (9) 

G R = (A— S fi ) _1 is, by definition, the retarded Green's 
function of the system under the influence of existence of 
leads. In Eq.(|5Jl we defined the effective incident wave 
for the boundary sites, g = b + BC _1 d whose nonzero 
elements are: 

s{n yp ,o-') = 2iSpp,8 <ja ,t pa sin^x^y) (10) 

from Ea.ljg^.we have cf> = CTM + C _1 DG fl g. By 
expanding this expression we obtain explicitly the scat- 
tering amplitudes: 
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= -<W<W"W + yV« E xL'(yp')G i? K p ,cr',^ p cr)[2^ )9S sin(fc^)xm(yp)] (11) 



III. GENERAL DISCUSSION OF RIGOROUS 
PROPERTY 

A. The Green's function formula 

Now let's calculate the multi-terminal transmission 
probability in the Landauer-Buttiker theory 0, T p ? a , = 

£ 

velocity of mode ml in lead p' . With Eg. ((TT|l . for 

P ^ p',we get: 



rpma 
p' m' a' 



V -^ L , where iy m / = ^2tysin(fc^/) is the 



•p' a' 



rpma 



m.m' 



Tr(r?G^ G R ,„) 



(12) 



where 



m ^ 

and G^ CT , and G A a , are spin-resolved retarded and ad- 
vanced Green's functions. This is just the most fre- 
quently used Green's function formula 1]. Here we have 
deduced it in a rigorous manner from the scattering wave 
description. We point out here that self energy T, R de- 
duced in many text books 1] have neglected contribution 
from evanescent modes, while our expression in Eq.JSJ) 
with longitudinal wave vector given by Eq. Q was rigor- 
ous. 



density of states arising from incident waves in lead p as 
following: 

ma 

4 E G R (y , a',n pyi a)T p (y 1 ,y 2 )G A (n py2 a,y'a') 



Vl,V2,cr 

i A p (y'a',y'a') 



(13) 



where we defined A p (y'a' ,y'a') as the partial/flowing 
spectral function arising from current incidence in lead 
P- 

Consider the equilibrium state when all reservoirs has 
the same chemical potential. Then, at any energy, the 
carriers is incident equally from all reservoirs and we 
should sum up the incident lead index p for scattering 
wave functions when calculating flowing density of states: 

p,m<r 

= ^ E G R (y'a\n yi a)T(y 1 ,y 2 )G A (n y2 a,y'a'jU) 



n vi > n V2 refer to all boundary sites in the system which 
has nearest hopping term to leads, and T — £r p . Fur- 
thermore, we can prove an important relation lj. Since 

[G R ] _1 - [G A ] _1 = - E fl = i£) I* = *r 



then. 



G A — G R = iG R TG A = iA 



B. Some rigorous property of equilibrium state 

1. flowing/partial and full spectral function 

Now, following Landauer's spirit 0, we assume the 
reservoirs connecting the leads at infinity will feed one- 
way moving particles to leads according to their own 
chemical potential. Let's normalize the scattering wave 
function ip pma so that there is one particle in the in- 
cident wave, i.e., we change e lk ™ a to e lkm ' r /yL, where 
L — > oo is the length of lead p. Meanwhile, the density of 
states in lead p is = tttt — ■ Now let's we add all 

incident channel (ma) of lead p with corresponding den- 
sity of states(DOS). Then, we get the partial(or flowing) 



where A = i(G R — G ) is the full spectral function which 
plays a role of a generalized density of states inside the 
SO system taking the presence of all leads into account. 
In an open system, it's not absolute clear how to define 
state and density of states. Here through our deduction 
we have made explicit the physical meaning of A from the 
wave function description. We emphasize the particular 
boundary condition used in the leads, i.e., there is a ther- 
mal reservoir with the same chemical potential /x, in the 
far end for every lead, in the Landauer-Buttiker sense. 
For non-equilibrium state, i.e., not all of the leads have 
same chemical potential and there will be net current in 
the leads and system. At that time,we are interested in 
flowing density A p rather than total density A. 

On the above we used only the diagonal elements of 
A p and A. It is convenient to define the non-diagonal 
elements as flowing spin density density matrix: 
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• pin 



AP(y'a, y"l3) =J2WP)A(y'a)l/Hv p 

ma 

= ^2 G R {y / a,ny lp a)T p ( yi ,y 2 )G A (n y2p cr,y"P) 

A( y 'a, v"p) = £ r s {y"(])i>s{y'a)i/hv pm 

p,m<T 

]T G R iy'a,n yip a)T( yi ,y 2 )G A (n V2p a,y"P) 
= i(G R (y'a,y"(3)-G A (y'a,y"[3)) (15) 



2. All T odd quantities vanishes in time reversal symmetric 

system 

Now let's prove that in the Landauer-Buttiker descrip- 
tion, for a time-reversal symmetric Hamiltonian, all T- 
odd quantities in equilibrium are zero: 

Consider a local quantity: 



(16) 



a/3 



Typical examples: (l)for spin density on site i, we 
choose O u = j5ij<T. (2)for normally defined charge cur- 
rent on link we set O 4 - 7 = — O ji = iel2 X 2 where 
^2x2 is the 2d unit matrix, (3)for normally defined spin 
current on link we have O y = — = In 
general we require O u = O^' for physical quantity. 

Under time reversal manipulation, T = ia y K, where 
K is the conjugate operator, O 4 - 7 will transform as O % ^o — > 



(— l) a+,3 0^*. There are two classes of physical quanti- 
ties, depending on their symmetry under time reversal 
operation, i.e., either T-even or T-odd quantities. Spin 
density and normally defined link currents are T-odd 
quantities. In contrast, the normally defined spin cur- 
rents are T-even quantity. 

Now, let's calculate the expectation value of O y in 
mesoscopic equilibrium state. By definition, this expec- 
tation value per energy interval is: 



prncr 



-Lj2lM3P,ia)0% + h.c] 

a/3 

^£[(G* - gA )^°% + ( GR ~ G A ) ia , ti O}U 

(17) 



2tt 



a/3 



where we have used the relation (G fl )f = G A . More- 
over,the matrices (G R ) and G A are related by time re- 
versal operation. Generally, by definition, we have: 

(TGgT-% J/3 = (-1) Q+/3 G«* {B) = Gf 0M {-B) 
Here, for generality of discussion, we have used magnetic 
field B to break T symmetry explicitly. Furthermore, 
let's suppose that the Hamiltonian is T invariant, then 



(TG T~ )ia,jp = {-l) a+ "G?l fp = Gf^ a (18) 

The four terms in Ea. (|17|l can be grouped into two 
pairs, 



a(3 a/3 



a/3 a/3 

r 



Now, suppose O is T odd, then, (— l) Q+/3 0^ = ~O l ^ . Thus we have 



y^}Gf ia Q^ g - G A a j pO j p a } +y^[ G A si ia O % i fi -I- G^jpO^] 

a/3 a/3 

a/3 a/3 

I 



evidently, each summation on the right hand side of Ea. (|19|) is real. So, the average value for physical quan- 
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tity in Eg. l|17[) should be zero. 

On the above we have proved that for any time rever- 
sal symmetric system, all T-odd quantities, such as spin 
density and local charge current, vanish in equilibrium 
state. This result seems to be a transparent truth. How- 
ever, up to our knowledge, though many people accept it 
to be true, no one has proved it rigorously in Landauer's 
mesoscopic structure. The corresponding fact for a closed 
system in equilibrium is trivial, because in closed system, 
each quantum state has its T partner state (except the T 
invariant states for which the average of T-odd operator 
should be zero), thus for a T-odd operator, the sum of 
average is zero. However, in our open system case, for 
equilibrium state, we have an ensemble of scattering wave 
functions incident from all leads at fermi energy while all 
the states are not manifestly T paired. 

Under finite terminal bias, the system is driven into 
non-equilibrium state. Due to above property, the value 
of T-odd quantity is determined by fermi surface prop- 
erty for small bias. This is in accordance with fermi liq- 
uid theory. But for T-even quantities, e.g., spin current, 
we cannot reach to such a conclusion. In fact, we need 
other symmetry considerations to interpret some T-even 
quantities to be fermi surface property 

For later use, let's consider a geometrically symmetric 
system with double terminals, as shown in FigEl Since 
the spin density vanishes in equilibrium state, according 
to Eq.{I7I)we have: 

< S a (x,y) > eq =< S a (x,y) >i + < S a (x,y) >_/= 

(20) 

where < S a (x,y) >i means partial/flowing average of 
spin density due to scattering states incident from lead 
1, and < S a {x,y) >_/ means flowing spin density due 
to scattering waves incident from lead 2. This relation- 
ship will be employed to derive an important symmetry 
properties in next section. 



p'th lead can be calculated as, 



lead 1 
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lead 2 




O X 
2DEG 







FIG. 2: schematic geometry of two-terminal structure 



3. equilibrium terminal spin current is always zero 

Next, for case of reference, we will give a explicit proof 
that terminal spin current also vanishes in equilibrium 
state(which is known by others irrespective of the T 
symmetry of the Hamiltonian. In linear transport regime, 
we neglect the energy dependence of scattering matrices, 
then the terminal spin current polarized in a direction in 



r = 

p 



T p% )V q ~ - T^,)V P ] 



+2(iC 



R%)V P } (21) 



Where we used symbol R to represent reflection prob- 
ability in lead p, i.e., we define Rf™ = etc for clarity. 
On the above equation, the summation is taken over all 
terminals except p. The second term in {• • • } describes 
contribution to spin current due to injection and reflec- 
tion wave in p'th lead. This term is zero for normal cur- 
rent, where summation of all spin indices is performed. 

For equilibrium state, we have V q =constant for all q. 
Let's firstly write down the following identities: 

ET pa , - R prT - R pa = —M 
L qo' ^pa ^pa lu p 

qa' 

I 

qa' 
I 

^ t Tp a + Rp a + Rp a = Alp 

qa' 
qa' 

(22) 

Here M p denotes the number of conducting mode in 
lead p. The first two equations are simple sum rules. 
They describe current conservation: one incident electron 
in lead p must go somewhere. The last two equations 
follows from the unitarity of scattering matrix. Or, we 
can deduce them from the first two since (B) = 

T p °,(-B),RPp% '(B) = RpZ,{-B) due to T transformation 
property of scattering amplitudes. Adding up all the 
above equations we get Ipf &q \ = 0. 

Thus, we have showed explicitly that terminal equi- 
librium spin currents are identically zero, irrespective of 
details of Hamiltonian of the system(even in the presence 
of magnetic field) and geometry structure. This result is 
parallel to the fact that terminal equilibrium charge cur- 
rent always vanish, while current density may be nonzero 
inside system when the system is placed under magnetic 
field. 

E.I.Rashba pointed out 8] that in a bulk 2DEG sys- 
tem with pure Rashba type of SO coupling, in-plane po- 
larized spin current may have non-zero value in equilib- 
rium. This simple result raised the problem of how to 
define spin current properly in SO systems. A much con- 
cerned problem is that whether such background spin 
current useful? From the practical point of view, the use- 
ful quantity is terminal spin current |3l| . In a conceived 
all-electric devise, we should add nonmagnetic contact to 
induct spin current out of SO system. However, accord- 
ing to the above rigorous result, the resulting equilibrium 
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spin current in added terminals around the SO system 
should always vanish and we cannot simply make use of 
it. How does spin current decay abruptly at the bound- 
ary between lead and bulk system is an interesting prob- 
lem to be clarified. Such an understanding may help to 
bridge the apparent gap between mesoscopic approaches 
and macroscopic approaches [24| (say, linear response ap- 
proach). 

When there is current flowing through the SO system, 
there will be induced spin current on attached leads. In 
the following we will discuss some symmetry property of 
such spin currents in a two terminal structure. We will 
reveal an emergent continuity property of spin currents 
due to symmetry of geometry and SO Hamiltonian. 

For later use, here we write down the following equa- 
tion for a two terminal system: 

where I^ieq) IS ^ ne equilibrium spin current in lead 2, 
and Zffi) an< ^ ^2(-j) denote spin current in lead 2 under 
the condition of fixed charge current I (flowing from lead 
1 — ► 2) and -I (flowing from lead 2 — » 1) respectively. 

IV. SOME SYMMETRY PROPERTY OF 
TYPICAL SPIN-ORBIT COUPLING SYSTEMS 

Now let's discuss the symmetry of typical Spin-orbit 
coupling system in which a SO bar is symmetrically at- 
tached by two ideal leads. Let's model the SO bar by a 
Rashba Hamiltonian, H so = a(a x k) ■ z, whose discrete 
version is just Eq.(J2J. We assume a rectangular geometry 
as shown in Fig 2. Firstly, let's analyze symmetry proper- 
ties of the model Hamiltonian. For pure system, it's easy 
to see from Fig. 2 and Eq.lQ), Eq.Q, that the Hamil- 
tonian has combined symmetry of real space reflection 
and spin space rotation: a y P x and o~ x P y , where P x and 
P y denotes real space reflection y =*> —y and x => — x, 
respectively, a x and a y acting on spinor wave function 
on every site are spin rotational manipulation around x 
axis and y axis, respectively. Moreover, when there is no 
external magnetic field, the Hamiltonian for the entire 
system is T symmetric. From these symmetries, we can 
obtain the following symmetry relations: 

<S x , z (x,y)>i = - < S x>z (x,-y) >i 

<S y {x,y)>i = < S y (x,-y) >i 

< S x (x,y) >/ = -<S x (-x,y)>i 

< S ViZ (x,y) >i = < S y , z (-x,y) >/ (24) 

< ■ • • >/ means taking ensemble average for given lon- 
gitudinal current / incident from lead 1. The first 
two lines are results of symmetry o~ y P x which has been 
known widely before [80l|. The second two, however, are 
results from a x P y and T symmetry together. From 
a x P y , we get:< S VlZ (x,y) >i= - < S v<z (-x,y) >_/,< 



S x (x,y) >/=< S x (—x,y) >-/. From time reversal sym- 
metry, we have Ea. (|20|l . Combining these two, we obtain 
the last two lines in Eq.J21}. 

Similar symmetry relations can be easily obtained for 
Dresselhaus SO coupling model[Tlj for which the discrete 
Hamiltonian reads: 

H° = t D ^m^ay* Ri+y -tf Ri a x * Ri+x )+h.c.} (25) 

Ri 

As can be easily seen, we have symmetry a x P x , o y P y 
and T in this Dresselhaus model. In this case, we get the 
following symmetry relations: 

< S y , z {x,y)>i = - < Sy, z (x,-y) >j 

<S x (x,y)>i = <S x {x,-y)>i 

<Sy(x,y)>! = -<S y (-x,y)>i 

<S x , z (x,y)>i = < S X)Z (-x,y) >i (26) 

In the presence of both Hf a and H® with geometry 
shown in Fig. 2, it can be checked easily that the total 
Hamiltonian is invariant under a z P x P y . Thus we have 
symmetry relations: 

<S z {x,y)>! = - < S z (-x, -y) >i 
<S XtV (x,y)>i = < S x , y (-x,-y) >i 

(27) 

Up to our knowledge, not all of the symmetry relations 
Ea.pijl. Eq.JSBJ and Eq. l(2Tjl have been reported previ- 
ously. In particular, Ea. (|27|l and the last two lines in 
Eq. (|24|l as well as in Eq. H26JI follow from T symmetry to- 
gether with geometric symmctry( spacial reflection com- 
bining a spin rotation manipulation) are firstly found in 
this paper. These relations provide important accuracy 
test for experimental measurements and numerical calcu- 
lations. For example, for a two-terminal structure with 
Rashba coupling in the middle SO bar, the theoretical 
assertion made in |l2j | that longitudinal spin component 
is insensitive to the reversal of the bias voltage is con- 
trary to our Ea. l2UI) . the numerical result 12] that spin 
density < S z (x,y) >v is not symmetric with x —* ~x is 
contrary to the last line of Ea. l|24(l . 

The issue of spin current has triggered hot discussion 
recently. The central problem lies in that spin is non- 
conservative in spin-orbit coupling system. In bulk sys- 
tems, proper definition and physical consequence of spin 
current are currently under debate. However, in a multi- 
terminal structure, spin current on the leads can be de- 
fined unambiguously. In the following let's analyze the 
symmetry of spin current in two terminal structure with 
geometric symmetry ( spacial reflection combining a spin 
rotation manipulation). The terminal spin current on a 
particular lead was a summation of link spin current (see 
discussion below Ea. 1)16(1 ^ in transverse direction. Then, 
for a system with both Hf a and H® present, by taking 
into account Eq.J33J) together with geometry symmetry 
and following similar argument in the derivation of sym- 
metry relations for spin density, we get a symmetry rela- 
tion similar to Eq. l|27|l : 
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dissipative systems. 



J 2(/) - 



•2(/) 



'!(/) 



(28) 



Such symmetry of terminal spin current is directly re- 
lated to symmetry of spin polarization Eq. I|27l) since the 
terminal spin current can be imagined to be injected from 
boundary spin density. The symmetry connection be- 
tween these two quantities is clearly depicted in Fig.J2J. 
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FIG. 3: schematic picture of symmetric relation between 
boundary spin density and terminal spin current in a two- 
terminal structure 

Interestingly, Ihn can be interpreted as conserved cur- 
rent because the quantity of it injecting into and leaving 
out of the SO region is identical. This emergent conti- 
nuity property, however, is the result of geometric sym- 
metry. In contrast to a recently proposed spin current 
dcfinition|23j| in bulk system, for which a coarse graining 
process(though physically beautiful, such process is hard 
to done rigorously) is necessary, the continuity of spin 
current in our case is due to dynamical and geometric 
symmetry. This interpretation has at least the advan- 
tage of physical transparency. The different property of 
terminal spin current polarized normal and parallel to 
the 2DEG plane is an interesting topic to discuss. As 
mentioned on the above, the equilibrium in-plane spin 
current is not useable in an all-electric device. Under 
the condition of flowing current in an ideal system, ter- 
minal in-plane spin current lacks a simple physical inter- 
pretation as a current while the out-of-plane spin current 
does. From these simple symmetry relations, we believe 
these components may play very different role in spin 
transportation in real systems. It is interesting to con- 
sider more geometries to further expose such symmetry 
differences and study the different role of in-plane and 
out-of-plane spin current to spin transportation in real 



FIG. 4: The charge and spin conductances as a function of 
Rashba coupling constant ta. The system is taken to be 100 x 
40 sized and Dresselhaus coupling constant is fixed to be 0.02t 

To have a quantitative impression on spin current, we 
performed some numerical calculation. Let's consider a 
two-terminal system with both Rashba and Dresselhaus 
coupling. We adopt a discrete representation and fix 
t jt = 0.02. The current is assumed to flow from lead 1 
to lead 2 for a system with 100 x 40 lattice sites and both 
leads are of the same width as system, the charge con- 
ductance a c and spin conductances a x ' v,z are shown in 
Fig.Q. The spin conductances are much smaller com- 
pared to charge conductances and both of them oscil- 
late rapidly in large tn/to regime. For tn/tr) < 2, a c 
decreases smoothly and a x ' v ' z are very small. When 
tn./tD > 2, both a c and a x s v ' z lines will show series of 
peak structures. The line-shape of a y s will follow closely 
to that of a c . On the other hand, the a x,z values will de- 
crease in the large tn regime due to symmetry property 
stated below. Once spin current become experimentally 
detectable, the peaks structures in the Fig.QJ are on the 
first place to be verified. 

For pure Rashba system, due to y — > — y symme- 
try(combining spin rotation) we have: 



jX.Z jx^z p. 

i 2(7) — — U 

TV _ TV 

J 2(I) - J 1(J) 

similarly, for pure Dresselhaus system, we have: 

jy,z _ jy,z _ „ 
1 2(I) — — U 

tx _ jx 

- J i(J) 



(29) 



(30) 



The symmetry relations Ea.p ^ .lj ^ .pC] )! can hold 
even in the presence of magnetic field (which keep geom- 
etry symmetry), since Ea. l|23|) is independent of T sym- 
metry. 
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On the above we have obtained some new symmetry 
relations relating to spin polarization as well as terminal 
spin currents for a two-terminal structure. Some other 
symmetries for a two-terminal waveguide in the pres- 
ence of magnetic field modulations is obtained in Ref . [32| . 
All these symmetry relations may find application in the 
study of spin transport theory. 

V. CONCLUSION 

To sum up, we have set up a detailed formulation 
to solve explicitly the scattering wave functions in a 
multi-terminal spin-orbit coupled system. We deduced 
some analytical properties in the scattering wave function 
description of mesoscopic physics following Landauer's 
spirit. In particular, we have (1) deduced rigorously the 



much used Green's function formula. (2) proved rigor- 
ously that all equilibrium value of T-odd quantities van- 
ish in the open multi-terminal structure. (3) reveal some 
new symmetry relations for two-terminal structure when 
the system has a Rashba or/and Dresselhaus form of SO 
coupling. These symmetry relations may provide accu- 
racy test to experimental measurements and symmetry 
restriction to numerical calculations. 
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